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DIRECT STEP METHOD FOR COMPUTING 
WATER-SURFACE PROFILES 


SYNOPSIS 

The direct step method proposed replace the computations trial and 
error used the standard step method for tracing water-surface profiles 
natural channels. semigraphical method Bernoulli’s theorem 
and the effect velocity-head changes. The effect eddy losses, 
bridge-pier losses, and overbank flow can also included when necessary. The 
proposed method, which conveniently applicable general, has been found 
the writer especially advantageous when several water-surface profiles 
must determined the same channel. 


INTRODUCTION 


Notation.—The letter symbols adopted for use this paper are defined 
where they first appear, the illustrations the text, and are arranged 
alphabetically for reference the Appendix. 

The proposed method seeks eliminate the trial-and-error computations 
necessary the standard step method for determining water-surface profiles 
natural channels. This accomplished indicating the hydraulic properties 
each section length river channel means graphs, such 
fashion that, knowing the water-surface level any section, the water-surface 
level the next section (upstream downstream) can obtained direct 


computation. 
THEORY 


represents the water-surface profile for varied flow between two sec- 
tions open channel. 

Because the total energy section (1) equal the total energy section 
(2) plus losses, 


comments are invited for publication; the last discussion should submitted 
September 1953. 


Structural Engr., Corps Engrs., San Francisco, Calif. 
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The losses the reach that are usually taken into account profile 
computations consist the loss head caused friction, the eddy losses 
caused enlargement contraction the channel, and bend losses. 

The friction loss the reach will computed Manning’s formula. 

The eddy losses (the part the velocity head that not recovered) are 
usually taken into account determinations water-surface profiles 


(a) They are included the value the rough- 
ness coefficient the Manning formula, and 
determined. This can done when the value 
has been determined from observations water- 
surface profiles the same similar reaches. The 
value this case will represent the total rate 
loss energy head including both eddy and friction 
losses. 

(b) The eddy losses are considered certain 
OPEN sections. The total loss head between two 
consecutive sections then will the sum the loss 

head caused the reach plus the eddy losses. 


Let the friction losses between sections (1) and (2) Fig. denoted 
Saver which the average the rates loss head due friction 
sections (1) and (2), and isthe length the reach measured along the center 
line between sections (1) and (2). 

Limiting consideration losses friction and eddy losses, the loss 
energy head between sections (1) and (2) Fig. 


the quantity parentheses the numerical difference between the 
larger and smaller velocity heads, irrespective which the larger, and and 
are the rates loss head due eddy losses are computed 
0.5 approximately. For gradually diverging reach may taken 
0.2 and for gradually converging reach, may taken from 0.1, 

When the effect eddy losses included the roughness coefficient 
Manning’s formula, and (as determined from Manning’s formula) repre- 
sent rates loss energy head sections (1) and (2). 

Substituting the foregoing expression for losses Eq. yields 


Corps Engineers Engineering Manual for Civil Works, Part CXIV, May, 1952, Chapter 
“Computation Backwater Curves River Paragraph 9-06. 
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For given discharge and given roughness coefficient 
formula yields 


1 


which the term always taken positive since the quantity 
parentheses merely the numerical difference explained previously. 
Thus, for constant discharge given roughness coefficient and given 


distance between sections, the expressions: and 
2 g A?, 2/2.2 A’, 


which the values and depend. The expressions the form: 


denoted and the expression the form: 


he 1 


when constant. 
DESCRIPTION PROCEDURE 


The direct step solution consists the graphical application For 
each chosen section the natural channel under consideration, the values 
and for range water-surface elevations are determined and tabulated. 
This procedure common both the standard step solution and the proposed 
direct step solution. least three water-surface elevations for each section 
are chosen, the range which the water-surface profile expected lie, and 
are used for computing the hydraulic properties. 


For given values and values and are determined each 


section for each these elevations, and the corresponding values and 


or 
(4) 


PROFILES 


“ 


< 


| 
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against for each section, shown Fig. For flows above critical depth, 
the water-surface profile determined upstream order. Starting with 
given water-surface elevation section, the value obtained 
from the appropriate curve. 

. V?; V2, . 

Neglecting the term when the effect eddy losses 
included the value the roughness coefficient can seen that the 
equal f(Z) for the next 23.0 
upstream section. Entering 
the curve for f(Z) for the 
next upstream section with 
this value, the corresponding 
water-surface elevation de- 
termined pro- 
cedure repeated from sec- 
tion section, tracing the 
desired water-surface profile. 
For flows below critical depth, 
the water-surface profile 
traced downstream and the 
procedure similar. Start- 
ing with the value f(Z) 
section and taking this 
for the next downstream sec- 
tion, the corresponding water- 
surface elevation deter- 

the eddy losses separately 
percentage the velocity- 
head change between the up- 18.0 
stream and downstream end 
sections reach, the pro- 
finding the water-surface ele- 
vations for two consecutive sections when proceeding upstream (as described pre- 
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are found for these water-surface elevations from previously plotted curves, 
shown Fig. The difference then added the value for 
the lower the two sections. This corrected value taken the f(Z) 
curve for the higher the two sections, and the corresponding corrected water- 
surface elevation determined. 


> 
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For flows below critical depth, the procedure similar, Eq. being appli- 

cable. However, the water-surface profile should traced downstream and the 

curve. 

Corrections for bend losses bridge-pier losses* reach may made 
manner similar that described for eddy losses. 

computing the values used plot the curve, the value 
downstream from that section should used. Similarly, computing 
for plotting the curve, the value upstream from that 
section should used. This insures that values will common sec- 
tions both ends any reach, that is, the same value used both 
sides Eq. choice the L-value important because sections taken 
natural channel usually are not equidistant from each other. 


correction for eddy losses 


TABLE Cross SECTIONS, FOR VALUES 


surface sectional Hygreule velocity, V, in feet 0.5 


(1) (2) 


Station 


(6) (8) 


0.00173 
9.57 A 0.00110 


wn 
~ 


+00 


+00 


+00 


POU 
CSC KOKO 


Com SON 


0.00102 


computing for any section, the effect overbank flow and 
varying can readily included the computations using weighted 
value the velocity for the sections when tabulating the properties 
for the purpose computing and after has been obtained 
for given value the corresponding values f(Z) and for any 

2 
other discharge can obtained multiplying the factor 
1 

The examples comprising the remainder this paper may help illustrate 

the use the proposed direct step method. 


No. 
discharge 33,500 per sec measured the San Lorenzo River 
Santa Cruz, Calif., let required determine the water-surface profile 


7 
238 
| 20.0 3500 248 
3970 258 
18.0 3580 265 
20.0 4100 273 
4630 282 
18.0 3080 248 
20.0 3550 258 
4050 268 
18.0 2650 225 
20.0 3100 235 
22.0 3550 245 
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station 18.5 above mean sea level for discharge 33,500 
persec. Assume that the value Manning’s 0.030, and that includes the 


effect eddy losses. 
Eq. applied graphically explained previously, neglecting the term 


The plotted cross sections from station station 


have been omitted from this presentation, but the hydraulic properties 


Station feet feet feet feet feet feet feet 
(5) (6) (7) (8) 


20.37 
21.70 
23.30 
19.58 
21.19 
22.92 
20.26 
21.67 
23.25 


eid 


for each section are tabulated Table for the range which the water-sur- 
face profile expected lie. 

The computations for and for three elevations are shown 
Table Curves f(Z) and were plotted against for each 
cross section, shown Fig. 

From Manning’s formula, was taken the length 
reach downstream from particular cross section; and was taken the 
length reach upstream from particular cross section. 

water-surface elevation 18.5 TABLE PROFILE, 
Fig. 20.70. Taking this 
the next upstream section sta- 
water-surface elevation found 20.70 
reached. The results the 
profile computations have been tabulated Table This method very 
advantageous the determination water-surface profiles natural 
channel, several profiles, each for different elevation the initial section, 


are required. 


+00 20.0 0.28 
22.0 0.19 
18.0 19.08 
+00 20.0 20.86 
22.0 22.69 
18.0 19.51 
+00 20.0 21,17 
22.0 22.91 
18.0 19.87 
+00 20.0 21.43 
22.0 23.13 
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this type problem are the tracing water-surface profiles 
tributary stream, for different stages the main river; river with 
tidal estuary, for different tidal elevations; channel connecting two 
reservoirs whose elevations vary. 


EXAMPLE No. 
Assume that required trace the water-surface profile from station 
per sec. The roughness coefficient Manning’s formula equal 


+00 


+00 


+00 


Si 


0.030, and the eddy losses included the value The water- 
surface elevation station 19.0 above mean sea level for value 

per sec from the corresponding values the previous problem, which 
33,500 per sec, multiplying them 1.43. The 

33,500 
values and obtained this manner have been tabulated Table 
and the corresponding values 
FOR VALUE been plotted Fig. 

Water-surface method. Fig. straight line 

(4) the f(Z) curve and the 

for station 00, the intercept 

between the line and 

the curve for f(Z) equal 

proportional dividers set the ratio 1.43:1, the intercept may set off 
the small end the dividers; and the intercept laid off with the large 
end shown—at the same water-surface elevation. The intercept there- 


‘ 
Station 
(1) (2) (3) (4) (6) 
22.43 
+00 20.0 2.43 22.43 
2.26 19.54 20.26 
20. 1.70 21.23 21.70 
1.32 22.99 23.32 
3.24 20.16 21.24 
20.0 2.39 21.67 22.39 
1.79 23.30 23.79 
22.0 
22.03 
20. 
23.62 
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fore will the intercept between the line and the curve for f(Z) for 
obtained similarly; and curves f(Z) and can obtained for 
each section for 40,000 per sec from the curves for 33,500 
per sec. The water-surface profile then may determined from this new set 
curves. The computations are tabulated Table 


No. 


For discharge 33,500 per sec the San Lorenzo River Santa 
Cruz, required determine the water-surface profile from station 
station For this discharge, the water-surface elevation station 
18.5 above mean sea level. The value Manning’s 0.030. 
Eddy losses are computed 50% the change velocity head when the 
velocity head decreasing downstream. Eddy losses may neglected when 
the velocity head increases downstream. 


AND 0.5 


losses 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 
80 +00 20.70 0.91 18.50 18.50 
85 +00 20.70 21.03 0.53 21.21 19.83 9.83 
89 +00 21.03 21.21 21.70 0.71 0.18 21.86 19.85 20.05 
94 +00 21.70 21.86 0.87 35 20.53 


All values are feet. 


The water-surface profile may determined application 
taking into account eddy losses. The value will 0.5, and the corre- 


2 
sponding terms are tabulated Table Fig. the curves for 


0.5 have been plotted and should used correct for eddy losses 
explained previously. The computations for determination the water- 
surface profile are tabulated Table Starting station with 
water-surface elevation 18.50 (see Col. Table 6), the value 
may found from the appropriate curve Fig. 20.70 (see Col. 


72 
Table 6). The corresponding value 0.91 (see Col. Table 6). 


Taking the value 20.70 the f(Z) curve for the next upstream section 


station (see Col. Table 6), the corresponding water-surface 
2 
elevation (shown Col. Table 19.83 ft, and the value 0.53 


(see Col. Table 6). the velocity increasing downstream and 
eddy losses are neglected, that the corrected water-surface elevation still 
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for water-surface elevation 19.83 found 21.03 ft. Going with this 
value the f(Z) curve for the next upstream section station 00, 
the water-surface elevation found 19.85 and the value 


2 
Therefore, the velocity head has decreased downstream from 


station station 00, and the eddy losses are 0.71 0.53 
0.18 ft. order include the effect these eddy losses, add 0.18 
21.03 which was previously found the value station 
00, yielding the value 21.21 ft. Going the f(Z) curve for station 
with this value, the water-surface elevation found 20.05 ft. 
This procedure repeated, accounting for eddy losses when necessary, 
until the water-surface profile has been computed over the required distance. 


APPENDIX. NOTATION 


The following symbols, adopted for use the paper and for the guidance 
discussers, conform essentially with ‘‘Letter Symbols for 
(ASA Z10.2—1942): 

area the cross section, square feet; 

acceleration due gravity, feet per second per second; 

ratio eddy losses velocity-head change between any two con- 

secutive sections; 

roughness coefficient the Manning formula; 

wetted perimeter, feet; 

discharge, cubic feet per second; 

hydraulic radius feet; 

rate loss head caused friction any section, determined 

the Manning formula; 

mean velocity flow any section, feet per second; and 

water-surface elevation any section, feet. 


